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Abstract 

The energy of the two-component Fermi gas with the s-wave contact interaction is 
a simple linear functional of its momentum distribution: 

internal = tfttC j \vam + ^(h 2 k 2 / 2m) (n ka - C/k A ) 

kcr 

where the external potential energy is not included, a is the scattering length, Q 
is the volume, n^ a is the average number of fermions with wave vector k and spin 
a, and C = lim^oo k 4 n^ = lim^oo k^n^i . This result is a universal identity. Its 
proof is facilitated by a novel mathematical idea, which might be of utility in dealing 
with ultraviolet divergences in quantum field theories. Other properties of this Fermi 
system, including the short-range structure of the one-body reduced density matrix 
and the pair correlation function, and the dimer-fermion scattering length, are also 
studied. 
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1 Introduction 



It is a textbook fact that the energy of a noninteracting two-component Fermi 
ffcLS IS cL simple functional of its momentum distribution: 



-^noninteracting ^ t ^k^kcrj 



where is the average number of fermions with wave vector k and spin a 
[Ij, 6k = /l 2 k 2 /2m, h is Plank's constant over 27r [2], and m is each fermion's 
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mass. This equation is valid for any nonrelativistic state; it is also valid when 
we have only a few fermions. 

What happens to the above simple formula if a contact interaction between 
the two components ("spin states") is turned on, characterized by a nonzero 
s-wave scattering length a? The system's wave function now behaves like 

= (l/r-l/a)A + 0(r) (1) 

if the distance r between two fermions in different spin states is small. Here A 
depends on the position of the center of mass of the two fermions, as well as the 
positions of the remaining fermions. This system has a superffuid phase with 
BEC-BCS crossover [3"fHlo] and other interesting properties, and has attracted 
a lot of experimental and theoretical work. So such a basic question deserves 
an answer. 

Because of the interaction, the momentum distribution now in general only 
decays like l//c 4 at large k [6], and the kinetic energy diverges [6]. In a real 
system, a contact interaction is impossible, and the divergence is eventually cut 
off at k ~ 1/tq, where r is the range of the interaction [7|. The interparticle 
interaction energy also depends sensitively on the physics at the short scale 
r . The sum of the two energies (denoted -External below), however, should be 
independent of such short-distance physics, when r is much smaller than the 
other relevant length scales in the problem, including \a\. 

We show in this paper that there is still a simple relation between the energy 
and the momentum distribution, which is independent of the details of the 
short-range interactions, except the scattering length a [8]: 

QC , k 2 ( C\ 

where 

C = lim k A n^ = lim k^n^, (3) 

k^oo k— >oo 

Vt is the volume of space j9], n ka = (cj^c^), and c k(T is the standard annihila- 
tion operator of a fermion with wave vector k and spin a. Equation ([2]) holds 
for any finite- energy states, no matter few-body or many-body, equilibrium 
or nonequilibrium, zero temperature or finite temperature, superfluid state or 
normal state. It holds for any populations of the two spin states (balanced or 
imbalanced). [ID] 

In the noninteracting limit, namely kpa — > C oc a 2 [IT] , and the first term 
on the right-hand side of ([2]) goes to 0. Here kp is the Fermi wave number. In 
the unitarity limit, in which kpa oo [12], C approaches some finite value, 
and the first term also vanishes. 
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Some technical remarks about Eq. 1) here and throughout this paper, 
"infinite momentum" physically stands for a momentum scale much lower 
than l/r , but much higher than any other relevant momentum scales in the 
problem, namely l/|a|, 1 /I, and I/A^b, where / is the typical interfermionic 
distance, and Xjb is the typical de Broglie wave length [13]; 2) -External is the 
expectation value of the total energy minus external potential energy |14j ; 3) 
the summation over momentum is understood as J2k = ^ / [d 3 &/(27r) 3 ], and 
there is no infrared divergence. 

In Sec. [2] a generalized function A(k), called the Lambda function, is intro- 
duced. Our motivation is to streamline the formulation of the s-wave contact 
interaction problem. Properties of A(k) are discussed. It is then used to study 
the two-body problem. 

In Sec. [3] another generalized function, L(k), is introduced. This and A(k) 
span a two dimensional linear space. A special element in this space, r/(k), is 
defined. 

In SecH]we prove ([2]) using r/(k). Physical implications are discussed, including 
a possible experimental test of ([2]), and results concerning the pair correlation 
strength and the one-body reduced density matrix. 

In Sec. we formulate the whole s-wave contact interaction problem (both 
few-body and many-body) in momentum space in a straightforward and un- 
ambiguous way, using the aforementioned generalized functions [15]. Using 
this formalism, we rederive ([2]) in a simpler way, and solve the three-body 
problem for the low energy scattering between a fermion and a weakly bound 
dimer of fermions. 

In Sec. [6] we summarize our findings, and discuss the possibility of extending 
the ideas contained in our method to some other important physical systems 
which appear at first sight to have little to do with the Fermi system considered 
in this paper. 



2 The Lambda Function 

2.1 Introducing A(k) 

If the interaction has a range tq that is much shorter than the other length 
scales in the problem (including \a\), it can be replaced by the famous pseu- 
dopotential [16J 

Ana d , A . 

V= 6(r)—r, 4 

m or 
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where r is the distance between two fermions in different spin states, and the 
partial derivative is carried out for fixed center-of-mass location of these two 
fermions, as well as locations of the remaining particles. When r — > but a 
is kept constant, the pseudopotential becomes exact. 

To formulate our current problem in momentum space, we perform a Fourier 
transform. To be concrete, we first consider two fermions interacting in vac- 
uum. The wave function associated with their relative motion, 0(r), satisfies 
a simple Schrodinger equation, 

E r 4> = -— 5(r)|-M, (5) 
m m or 

where E r is the energy of relative motion. Writing 0(r) = / ^p-0(k)e lk ' r , we 
get 



/ Wf{ Er ~ ^ (kyk ' r = ^ 1 ^^{^(r^). (6) 



Ana f d 3 k ~ d 



5(r)^-(re ikr ) = 5(r)A(k), (7) 



To proceed, we define 



dr 

where A(k) is so far an unknown function of k. 
Clearly 

A(k) = 1 (ifjfe<oo), (8a) 

but also, if we multiply both sides of Eq. ([7j) by l/k 2 , and integrate over k, 
the left side becomes oc <5(r)J^(r/r) = 0; we are thus led to an equality 

d 3 fc^ = 0, (8b) 

in apparent contradiction with Eq. (!8a|) . This difficulty has hampered the 
direct application of the pseudopotential in momentum space. 

Actually there is no contradiction here at all. 

Contradiction can only arise if we force a third equality, namely the integral 
of A(k)/k 2 over the whole k-space is equal to the limit of the integral over a 
finite k-space region which expands without bound. If we insist on both (IHaj) 
and (I8bl) . as we should, we must give up this third equality. 



At first sight, this decision is alarming, since the full-space integrals of all 
known functions are defined in terms of such a limit. In fact this is the standard 
definition of full-space integrals in mathematical textbooks. However, if A(k) 



4 



is accepted as a special generalized function, A(k)//c 2 does not have to obey 
this rule. 



For convenience we postulate one more property for A(k): 

A(-k)=A(k), 
whose usefulness will be clear shortly. 



Equations (l8a| . (l8bl and (IScjl define a generalized function, A(k), where k is 
a three-dimensional vector. 

Equation (j7j) is then a corollary of such a definition. 



Mathematical properties of A(k) 



Because A(ck) (where c is any nonzero finite real constant) satisfies the same 
three basic equations as A(k) itself, 



A(ck) = A(k). 

Similarly 

A*(k) = A(k). 
For any finite constant vector ko, 



d 3 A; 



A(k) 



(k - k ) 



(9a) 



(9b) 



(9c) 



To prove (l9cl) . we rewrite its left hand side as 



d 3 fcA(k) 



1 



(k 



ko) 2 



+ 



A- 2 



by using (I8bj) . Here c is any finite constant. If c = — 1, the integrand decays 
like 1/k 3 at large k, in a given direction. Without Eq. flHcl) . we could not 
determine the integral unambiguously. Using ([8c]) . however, we can take the 
average between the expression in the bracket and its spatial inversion, and 
get 

1 1 c 
1 1 

k 2 



d 3 fcA(k) 



2(k-k ) 2 



2(-k-k ) 2 



If c = — 1, the integrand now decays like l/k A at large k, and the integral is 
dominated by a finite region of the k-space (whose size is of the order fc ), so 
according to Eq. (l8~aj) . the Lambda function can now be dropped. The resultant 
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ordinary integral turns out to be 0. Equation (|9c|) is thus proved. It follows 
that for any finite constant k and nonzero finite real constant c 



A(ck-ko) =A(k). (9d) 
Similarly, it is easy to show, eg, 

d 3 k A(k) a 



(2vr) 3 (k-k ) 2 + a 2 4tt 



(a > 0), (9e) 



which is the k-space representation of the simple fact that the symmetric 
average of e~ ar+ °' r /47rr at r — > 0, excluding the l/Anr term, is — a/An. 
Equation (l9"e|) is a generalization of (I9c|) . 

A notable corollary of the Lambda function, as illustrated by fl9e|) . is that the 
full-space integral of a function that is positive for all finite k's may still be 
negative. We will see that this phenomenon is actually useful (in Sec. I4.3p . 

The integrals involving the Lambda function still satisfy many familiar rules. 
The region of integration, for example, can be freely divided into some subre- 
gions (with one restriction; see below), and the total integral equals the sum 
of the integrals over these subregions. For example, 

k 2 Jk- 



k<K K" Jk>K k 2 

= (+AttK) + (-AttK) = 0, 

for any finite K > 0. In the ball region, k is finite and A(k) = 1, while in the 
external region the integral is identical with that of A(k)6(k — K) / k 2 over the 
whole k-space, which can be computed with the same method as is used to 
derive Eq. fl9cl) . Here for simplicity we separate the subregions with a sphere, 
but other arbitrary shapes of the boundaries are equally permissible, provided 
that infinite momentum is contained by only one subregion. 

Like the delta function, A(k) can freely participate in operations like addi- 
tion, multiplication, integration, and Fourier transformation. It also commutes 
with other functions and quantum mechanical state vectors and operators, ie, 
A(k)X = XA(k), where X may or may not depend on k. 



2.3 The Fourier transform o/A(k) 



The Fourier transform of A(k) is 

d 3 k 



A(ri " / (i^ A(k)eik ' r - (10) 
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Since A(k) is both even and real, so is A(r): 

A(-r) = A(r), (11a) 
A*(r) = A(r). (lib) 

It follows from Eqs. (JSaJ) and (HQ} that 

A(cr) = |cr 3 A(r), (11c) 

a property reminiscent of the three-dimensional delta function. Here c is any 
nonzero finite real constant. 

Equation (19c]) can be expressed in r-space as / d 3 rA(r)e lk °' r /r = 0. Expanding 
this result in powers of k , we get 

/d 3 rA(r)- = 0, (lid) 
J r 

Jd 3 rX(r)v = 0. (lie) 

The fact that A(k) is equal to 1 for any finite k leads to 

J d 3 rA(r)/(r) = /(0), (llf) 

for any function / which is finite and continuous in a region containing r = 0. 
This entails, in particular, that 

A(r) = (ifr^O). (llg) 
From Eqs. (fTTdl) . (fTTel) . and (|TTfj) . we get 



/d 3 rA(r)^ = 0, (llh) 
J r 

for any function g which is smooth in a region containing r = 0. This equation 
becomes obvious if we write g(r) = A + B • r + 0(r 2 ). 



According to Eq. (llg), all the integrals involving A(r), like the ones showed 



above, can be restricted to a neighborhood of the origin without affecting 
their values. This is an important similarity between A(r) and 5(r). The two 
generalized functions mainly differ in two aspects: 1) J d 3 rA(r)/r = but 
/ d 3 r<5(r)/r is divergent; 2) / d 3 rA(r)f = but / d 3 r<5(r)r is undefined. 

A(r) is related to the operator 5(r)(d/dr)r: 





J d 3 r5(r) — [r/(r)] = J d 3 rA(r)/(r), (Hi) 
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for any function / for which / d 3 r<5(r) J^[r/(r)] is well-defined. 

Finally, if we integrate both sides of Eq. (J7J) over r, and use Eq. flllil) . we get 

J d 3 rA(r)e ik ' r = A(k), (llj) 

which is just the Fourier transformation from the r-space to the k-space. This 
equation is consistent with (fTUI) . because the Lambda function is even. 

It will be clear that A(k) and A(r) allow for a much more flexible treatment 
of the s-wave contact interaction problem than the operator 5(r)(d/dr)r. 



2.4 Two-body problem 

The right-hand side of (jSJ) equals / ■^4 3 'A(k / )0(k')5(r) according to (0). Ex- 
panding 8(r) in terms of plane waves, and comparing both sides of Eq. ([6]), 
we obtain the k-space representation of (jSJ), 

E r m = k -m + ^ / |^A(k')0(ko (12) 

or, equivalently, 

w) = ^»+^i;AM«k'). 

To familiarize ourselves with this formalism, we consider a simple exact solu- 
tion to (fl2"j) . namely the bound state (E r < 0). Let / denote the second term 
on the right-hand side of ( fl~2l) . Solving (112j) formally, we get 

0(k) = -f/(k 2 /m - Er) (f^O) (13) 

which is smooth for all k. Substituting ( fTBl back into the definition of /, and 
using (l9ell . we get / = fa^—mE r . So a > 0, E r = —1/ma 2 , and 0(k) oc 
l/(k 2 + a -2 ), in perfect agreement with the established wisdom |17j . 

We can easily extend the above approach to three or more particles (see 
Sec. 15.21 for an illustration). 
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3 Short-Range Selectors 



3.1 L(k) 



If the connection between / d 3 k and hin^^oo I\k\<K d 3 ^ is not universal, noth- 
ing can prevent the existence of another generalized function, L(k), defined 
as follows: 

L(k) = (if k < oo), (14a) 
d 3 k L(k) . . 

= l - (14b) 

L(-k) = L(k). (14c) 

As generalized functions, A(k) and L(k) can both be approached by ordinary 
functions. If k is a wave vector, L(k)'s dimension is length, while A(k) is 
dimensionless. 



3.2 Properties of the L function 



I p/W^)^^) ( 15a ) 

for any ordinary function /(k). 

We now list some other properties of L(k); their proofs are similar to those in 
Sec. El 

L(ck) = |c| _1 L(k) (15b) 

for any real constant c ^ 0. 

L*(k) = L(k). (15c) 

L(k - k ) = L(k) (15d) 

for any constant vector ko. An integral involving L(k), like that involving A(k), 
can be freely divided into many subintegrals (provided that infinite momentum 
is contained by only one subregion). L(k) can also freely participate in various 
operations, and it commutes with other objects. 

Now turn to the coordinate representation of the L function, 
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For any ordinary function /(r), 

/ d 3 r/(r)/(r) = lim 4?rr/(r). (17a) 

J r^O 

Other properties of l(r) are listed below. 

Z(r) = (if r^O). (17b) 

i*(r) = l(-r) = l(r). (17c) 
/( cr ) = |c|- 2 /(r) (17d) 
for any real constant c ^ 0. Equation (|17dl) is consistent with the dimension 



of Z(r), namely [length] 2 . 

"d 3 r/(r)e- ik - r = L(k). (17e) 



5.5 Introducing the short-range selectors 



Let /i(k) be any ordinary function satisfying / 7§^-/i(k) = 1. Let /^(k) 
1/fc 2 + r'(k) be any ordinary function satisfying / ^p-r'(k) = 0. Let Si(k) 
A(k) and s 2 (k) = L(k). We have 



/ 7^3 S K k )/i( k ) = V ( 18 ) 



/i(k) and / 2 (k) span a two-dimensional linear space T. A(k) and L(k) span 
another two-dimensional linear space 5. Equation (fT8|) states that and 5 
are dual linear spaces. 

We shall call the elements of S (short-range) selectors, because for any function 
/(k) = J2i c i/j(k), we can selectively extract the coefficient q using an element 
in 5:/^ a :(k)/(k) = c i . 

All these functions can be Fourier-transformed to the r-space. Let §j(r) = 
/ (g^(k)e ik - r , /,(r) = / 0/,(k)e ik - r . 

'd 3 rs*(r)£(r) = <%, (19) 

which represents a group of nontrivial equations: 

Jd 3 rA(r) = l, Jd 3 rA(r)/(47rr) = 0, 
Jd 3 r/(r) = 0, Jd 3 r/(r)/(47rr) = 1. 
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Also, 



s,(r)=0 (ifr^O). 



(21) 



Each short-range selector corresponds to a linear functional which extracts a 
short-range property of an ordinary function /(r). 

At this point, we have completed a tool which will free us from ill-defined 
ultraviolet divergences, for the two-component Fermi gas with s-wave con- 
tact interaction [I]. In this formalism, no ad hoc large momentum cut-offs or 
dimensional regularizations are needed. 

Here we have presented our approach in a generic form, so that it might be 
possible to extend it to other physical systems involving contact interactions. 
For the other systems we may need more than two linearly independent short- 
range selectors. For instance, we may need three independent selectors when 
the two components of the Fermi gas with s-wave contact interaction have a 
mass ratio exceeding 13.6, so that there is Efimov effect which introduces an 
additional parameter for the interaction [HUH]. 



3.4 The rj-selector 



The ^-selector is a particular element in S: 

,,(k) = A(k) + ^ , (22a) 

This selector will play a crucial role in the s-wave contact interaction problem, 
because it selectively annihilates the relative wave function of two particles 
with such interaction: 



j d 3 r^(r) [1/r - 1/a + 0(r)] = 0. (23a) 



On the other hand, 

J d 3 r77(r)/(r) = /(0), (23b) 

for any ordinary function /(r) that is continuous in the neighborhood of the 
origin, so fj(r) behaves like the delta function for nonsingular functions. More 
properties of the 77-selector are listed below. 

j d 3 rr/(r)f = 0. (23c) 

77(r) = (ifr^O). (23d) 
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Tj*(r) = f](—r) = T](r). 
77 (k) = 1 (if k < oo). 
d 3 A; 77(k) 1 



d 3 k 
(2tt) j 



^(k)/(k) 



(2tt) 3 k 2 
+ lim 



AtTO K-*oo J\\s\<K (27r) 3 



47TCZ 

d 3 A; 



where /(k) is an ordinary function and c = lim^oo k 2 f(\t). 



(23e) 
(23f) 

(23g) 
(23h) 



4 Energy Theorem 



^ . 1 Mathematical formulation 



Theorem. If the system of fermions of equal mass m populating two spin 
states with s-wave contact interaction and scattering length a is in a smooth 
external potential V ext (r), and is in a state p = c^il^i) (where {4>i\4>j) = 
5ij, act > 0, and J2i — 1) satisfying two conditions: firstly each \<f>i) is a linear 
combination of energy eigenstates with coefficients of the combination decay- 
ing sufficiently fast at large energy such that the wave function 0j in the coor- 
dinate representation has no singularities other than those introduced by the 
interfermionic interaction, and secondly the probability aij decays sufficiently 
fast at large i such that 

C = X)a i C < , (24) 

i 

where C [defined by fl3])] is associated with the state p, and associated with 
\<f>i) |2D], then the system's energy expectation value is 

E = J2v(^n^ + J2f d 3 r\/ cxt (r) P(T (r). (25) 

Here r/(k) is defined in the previous section, is the momentum distribution, 
and Pair) is the spatial density distribution. 

Proof. The second term on the right-hand side of fl25l) is trivial. We will concen- 
trate on the first term, -External, which is physically the sum of the total kinetic 
energy and the interfermionic interaction energy, both of which are divergent 
in the zero-range interaction limit. However, -External can be unambiguously 
determined in this limit (with a^O fixed), using the pseudopotential method. 

Let us first consider the case in which the system is in a pure state |0) having 
exactly N fermions in the spin up state and M fermions in the spin down 
state: 
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NlW. I d3ri ' ' ' d3rArd3,Sl ' ' ' d3s Af0(ri, • • • , rjv, si, • • • , s M ) 

x V4(r x ) • • ^{(r^)^l(si) • • ^i(sAf)lO), (26) 



where |0) is the particle vacuum, and ^ CT (r) is the standard fermion anni- 
hilation operator at spin state a and spatial location r. The wave function 
4>(r 1 , • ■ ■ , r N , si, - • • , s M ) is completely antisymmetric under the exchange of 
any two fermions in the same spin state; it has also been properly normalized: 

y'd 3 r 1 ---d 3 r J vd 3 s 1 ---d 3 SAf|0(r 1 ,...,r J v,s 1 ,...,SA/)| 2 = 1, (27) 
so that = 1. 

The only role of the pseudopotential is to exactly cancel the delta function sin- 
gularities arising when the kinetic energy operator (which is just a 3N + 3M- 
dimensional Laplace operator divided by —2m) acts on the wave function. 
These delta functions arise when two fermions with opposite spins come to- 
gether. To understand this, we may examine Eq. (jSJ) closely; note also that 
this cancellation mechanism is carried over to the arbitrary-body cases, with 
or without external potentials. The internal energy of the fermions is therefore 
given by a simple expression, 

- 2m£ internal = -i- lim / d m+3M R 0*(R) V 2 0(R), (28) 
JS\M\ £^oJp( e ) 

where R is the shorthand for the 3iV + 3M coordinates of the fermions, and 
V 2 is the 3iV + 3M-dimensional Laplace operator. T>(e) is a subset of the 
3iV + 3M-dimensional configuration space, excluding the regions in which any 
two fermions with opposite spins have a distance less than e. 

Now define another quantity X, 

- 2mX = J d 3N+3M RdH rKt)0*(n ■ • ■ r N s, ■ ■ ■ s M ) 

r N 

0( r i • • • r i + t, r i+1 ■ ■ ■ r N s 1 ■ ■ ■ s M ) 



i=l 

M 

+ X! ^( ri ' ' ' r N s l " ' ' s i-l) s j + t, Sj + i ■ ■ ■ Sm) 
3=1 



(29) 



and we want to prove that X = internal- To do so, we divide the 3iV + 3M- 
dimensional R-space in this new integral into T>(e) and X(e), where T>(e) is the 
same as above, and X(e) is complementary to T>(e). Clearly, in the subregion 
T>(e), the integral over R is finite and continuous, and the result is a continuous 
function of t for t < e; then, according to the basic property of ^(t), it can be 
treated as the delta function, and we immediately see the integral is exactly 
equal to the integral in internal- 
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We then only need to show that the integral in X(e) approaches zero as e — > 0. 
Since the volume of X(e) is proportional to e 3 when e is sufficiently small, 
X(e) can be further divided into NM subregions, in one of which |r^ — Si| < 
e. All the other NM — 1 subregions are the same as this, due to fermionic 
symmetry We have omitted subregions with volumes of higher orders in e. 
In the thermodynamic limit, however, even when e is very small, we can still 
have many pairs of fermions, and the distance between two fermions in each 
of these pairs is smaller than e. But we note that these pairs are far apart if e 
is small, so they can be treated independently. 

So now we only consider the case in which |r x — si| < e, which is representative 
of the general situation. In this case, only two terms in the big bracket on the 
right side of Eq. (|29|) have the possibility of making contributions to the total 
integral which do not approach zero. One of them is the term in which ri is 
replaced by r x + 1, the other being the one in which si is replaced by S! + t. 
We only discuss the first term, since the logic is the same for the second one. 

To treat this term we make a coordinate transformation: r = r x — Si and 
r = (r x +Si)/2, and represent r 2 , • • • , r N , s 2 , • • • , s M with a single 3N+3M—6 
dimensional vector R'. We then first do integral over r (and r < e), then do 
the integral over t, and finally integrate over ro and R'. 

Expanding 0(R', r , r) in this A(R', r )(l/r - l/a) + 0(r) (according 

to the short range boundary condition), we write our target integral in the 
form 

Y = J d m+3M - 6 R'd 3 r [ d 3 ^(t)V t 2 ^(R',r ,t), 

where we have omitted the constant coefficient since it is irrelevant to our 
question, and 

#(R',r ,t) = f d 3 r[A*(R',r )(l/r-l/a) + O(r) 

Jr<e 1 

x ~A(R', r + t/2)(l/|r + t| - l/a) + 0(|r + t|) 



K should be expanded in powers of the small t, before we can carry out 
the integral over t; t is regarded as much smaller than e, because rj(t) is 
zero for any nonzero t. In such an expansion, any term which contains a 
factor e raised to any positive power should be omitted. Also, any term which 
is of the order t 3 or higher should be omitted, since it contributes nothing 
to the integral /d 3 t^(t)V?[-]. We then have K ~ c x (t/2 - t 2 /Qa) + c 2 ■ it 2 , 
where a = -4vr|A(R, r )| 2 , and c 2 = -irA*(R', r ) V ro A(R, r Q ). But now 
V\K ~ c\{l/t — l/a) + 4c 2 • t, and according to the short range properties of 
rj(t), the integral over t vanishes. 

The above analysis shows that the contribution to X from X(e) approaches 
zero as e -> 0, so X = E inteTnal . 
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Equation (129|) can be easily rewritten in the second quantized form: 
-2mX = / d 3 rd 3 ^(t)V t 2 ^(r)^(r + t)|0). 

Expanding i^ a {v) = Vt~ l l 2 X) kl c kio- exp(iki • r) and similarly for ^(r + 1), and 
carrying out the integration over r and t, we get 

-2mX = MErtk)(-* 2 )4<*>>- 

kcr 

So £ internal = Ek^(k)(A; 2 /2m)(0|4 (T c k(J |0). 

So far, we have proved (125]) for a pure quantum state, with fixed numbers of 
fermions in the two spin states. 

If the pure state is not an eigenstate of particle numbers in the two spin states, 
we can expand it as a superposition of such eigenstates; since the interaction 
conserves the number of fermions in each spin state, the off-diagonal matrix 
elements of the Hamiltonian between these eigenstates are zero, and the the- 
orem remains valid. 

If the system is in a mixed state, described by the density operator p, and the 
different eigenstates of p independently satisfy the theorem, then the statistical 
ensemble of these states still satisfies the theorem, provided that (pHl) holds. 

Using Eq. (pEj) to reexpress the first term on the right-hand side of ( 1251) . we 
get ©. 

From the many-body wave function, we can also prove Eq. ([3]) by expand- 
ing (ipl(r)if) a (r + t)) at small t. In this expansion there is a singular term 
(proportional to t but independent of t) that is independent of a: 

W*(r)^(r + t)> = p„(r) - ^ + ip„(r) • t + 0(( 2 ), (30) 

C « " { N-mM-iy. K +3M - 6 *'»)| 2 , (3D 

where A(R!, r ) is defined in the above proof. Consequently nk CT decays like 
C/k A at large k, and 

QC = J C(r)d 3 r = J. (32) 
We can also show a result for the pair correlations: 

(p T (r - t/2)p i (r + t/2)) = ^(i - £) + 0(f) (33) 

at small t, where p CT (r) = ipl(r)ip a (r). 
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4-2 Other proofs 



For those who are not satisfied with the proof of the energy theorem with 
the 77-selector, we have in principle at least two other proofs. One of them 
is almost identical with the proof presented; instead of directly using the 77- 
selector to annihilate some terms, we may do a detailed analysis of the behavior 
of ( , 0^(r)'0 o .(r + t)) at small t, and, besides the term proportional to t (inde- 
pendent of the direction of t), we also need to analyze the term proportional 
to t 2 (also independent of t). There is also a term proportional to the inner 
product of t 2 t and a direction vector, but it contributes nothing if we do the 
symmetric integral in momentum space (as is done in Eq. (T5])). 

The other proof is inspired by what Eq. (T5]) tells us. Instead of studying an 
idealized s-wave contact interaction model, we may study a short-range (ro) 
attractive interaction potential, and fine-tune the depth of the potential to 
achieve a specified scattering length (|a| ^> r ). Then, we can divide the total 
internal energy in two pieces. One of them is the integral of kinetic energy 
up to a momentum scale K, where K is much smaller than l/r but much 
larger than the other characteristic momentum scales in the problem. The 
other piece is the kinetic energy integral from K to much higher than 1/Vo, 
plus all the interfermionic interaction energy. It is then possible to show that 
the second piece can be expressed in terms of C ~ i^ 4 n K(7 : approximately 
[n/(2a) — K]QC/(27r 2 m) + 0(1/K), with other errors that vanish in the limit 
r — > (but the scattering length is kept constant). In the two-body case, 
this approach can work out without much difficulty In the many-body cases, 
it is however tricky to give a rigorous proof; but there is a heuristic physical 
picture: at large K, we are effectively probing those fermions each of which is 
close to another fermion (with distance ~ 1/K), and for each of such fermion 
pair, two-body physics is a good approximation. 



4-3 Physical Implications 

One can measure the momentum distribution experimentally. A well- 
known method is to suddenly switch off both the interaction between fermions 
and the external confinement potential, to allow the fermionic cloud to expand 
ballistically. Because particles with different momenta move at different veloc- 
ities, eventually the spatial density distribution will reflect the momentum 
distribution. (The spatial distribution can be measured with some imaging 
technique.) See, eg, Ref. [ZT] . 

From the measured rZk CT , one can compute a partial kinetic energy, by summing 
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Fig. 1. Partial kinetic energy T(K) versus K. The vertical intercept of the asymptote 
at Kq = ir /2a equals the internal energy, according to ([35]) . 

contributions from all k's up to scale K: 

wo = £ ,34) 

|k|<K,<7 

and plot T(K) versus K. Because of the nonzero C [22] (in the C/k A tail 
of the momentum distribution), T(K) does not approach any finite value for 
K <C l/ r 0) but instead, it approaches a straight line with a positive slope, 
which is just the asymptote of the curve T(K). What Eq. (J2J) states is simply 
that the internal energy is exactly equal to the vertical intercept of such an 
asymptote at K = 7r/(2a). See Fig. [1] for illustration. More specifically, 

T(K) = £ internal + (K-pi+ 0(1/K), (35) 

2-K z m \ la J 

for K much larger than the other momentum scales in the problem but much 
smaller than 1/tq. This result is universally valid, as is stressed in Sec. [TJ 

The 0(1 /K) term on the right-hand side of (|35|) is related to the second 
order singularity of the relative wave function of two fermions with opposite 
spins. The first order singularity is like 1/r and the second order is like r. The 
interference between these two orders gives rise to a term oc 1/K at large K. 

If we plot T(K) all the way up to the scale 1/ro, the above asymptote behavior 
breaks down, and the curve eventually approaches a horizontal line, associated 
with the total kinetic energy, which is much larger than the internal energy. 

It is simple mathematics to confirm the energy theorem in the case of an 
isolated bound state of two fermions, whose momentum distribution is [17J 

nka (|k-k | 2 + a- 2 ) 2 ' an' 
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where 2ftko is the total momentum. In this case the first term on the right-hand 
side of Eq. (j2J) equals +2h 2 / (ma 2 ), and the second term equals —3h 2 / (ma 2 ) + 
(hk ) 2 /m. The total energy is —h 2 /(ma 2 ) + (hk ) 2 /m, simply a sum of the 
binding energy [T7] and the center-of-mass kinetic energy. Note that the total 
energy is negative if ko < a -1 . 

Here we see that the symbol r/(k) has to break the conventional law of integral 
(that the integral over the whole k-space equals the limit of integrals over 
finite regions of this space), otherwise it would be impossible for the first term 
on the right-hand side of Eq. (|25[) to generate any negative number, since 
77 (k) is equal to +1 for any finite k and the momentum distribution is also 
everywhere positive. The negative binding energy is associated with this purely 
positive momentum distribution. We could of course reject this modified notion 
of integration, and stick to Eq. (J2J) only, but then the beautiful structure of 
the problem, as shown by (I25j) . would be obscured. 

The energy theorem takes a particularly simple form in the unitarity limit: 
when a — > 00, the internal energy is just the intercept of the aforementioned 
asymptote at zero momentum. 

The quantum few-body problem in the unitarity limit in a harmonic trap is 
preliminarily studied in [23]; we may also study the momentum distributions 
of these few-body systems and relate them to the energies, using Eq. (135]) . 

Equations fl32l) and (133]) imply that the expectation of the number of pairs of 
fermions with diameters smaller than a small distance s is 

A^pair = ^ (36) 

if s is still much larger than tq. More specifically, the expectation of the number 
of pairs of fermions located in a volume element d 3 r, with diameters smaller 
than a small distance s is 

diV pair = ^d 3 r. (37) 

For these reasons, we shall call C(r) local contact intensity, X — J C(r)d 3 r = 
flC integrated contact intensity, and C average contact intensity (over volume 
tt). 

In a trapped Fermi gas, when the Thomas-Fermi approximation is valid, we 
may also have an approximate concept of local momentum distribution, valid 
within the limit set by Heisenberg uncertainty principle. Here we are using 
the term "Thomas-Fermi approximation" in its broadest sense, namely the 
fermionic cloud is divided into many portions: each portion is both so small 
that it is roughly uniform, and so large that it can be treated as in the ther- 
modynamic limit, so that all the thermodynamic quantities are meaningful. 
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Even when the Thomas-Fermi approximation is invalid (eg in few-body sys- 
tems), the local contact intensity C(r) is still an exact concept. The resolution 
8k of the momentum distribution within a local spatial region is of the order 
the inverse size of the region. In the large momentum part, where we can use 
a large unit for momentum, 5k is negligible, and a tail C(r)/k 4 is well defined. 

We also have a concept of local internal energy density eternal, inspired by an 
intermediate step in the proof of the energy theorem: 

6int ema i(r) = -E/d 3 * f(B)^(#Wr + s)) ( (38) 
and the total internal energy is exactly 

-^internal = / einternal(r)d 3 r. (39) 



The above results concerning local pair correlations and local energy densities 
can help us to address a fundamental problem for a finite number of fermions, 
namely how to construct a systematically improvable perturbation theory, 
with the Thomas-Fermi approximation as the zeroth order approximation. 
The problem is important for two reasons: 1) ultracold atomic Fermi gases 
realized so far typically contain hundreds of thousands of atoms, for which the 
Thomas-Fermi approximation has detectable errors, and 2) if we want to use 
a quantum few-body system to simulate the many-body thermodynamic limit 
(see, eg, [23J ) , finite size corrections must be taken into account. In fact few- 
body calculations are all that we can do. BCS theory, for example, is all about 
the two-fermion correlations in the presence of a many-body mean field. Quan- 
tum Monte-Carlo simulations are restricted to a small number of particles. In 
fact we can not solve Schrodinger equation for more than a handful of particles 
(except for the noninteracting systems or some very special systems). 

In the thermodynamic limit, the behavior of C or the complete momen- 
tum distribution on the entire (— 1/fcpa, T /Tp) plane (for balanced popula- 
tions of the two spin states) is worth studying. Here T is temperature, and 
Tp = h 2 k F /2mkB the Fermi temperature. Since the energy does not change 
smoothly across the phase transition line, the energy theorem implies that the 
momentum distribution also has some unsmooth change. 

Another implication of the energy theorem concerns a common dynamic pro- 
cess: the Fermi gas is initially confined in a trap, and then the trap potential 
is suddenly turned off (but the scattering length is not changed), so that the 
gas expands in the presence of interactions. During the expansion, the contact 
intensity continuously decreases, in such a way that the vertical intercept [at 
Kq = 7r/ (2a)] of the asymptote of the function T{K) remains constant because 
of energy conservation. This is a constraint on the intantaneous momentum 
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distributions. 



If the Fermi gas is in a motion (or if we are looking at a local portion of the 
expanding cloud in the local approximation), and the momentum distribution 
is roughly symmetric around k = k , it is better to use the following alternative 
formula to evaluate the internal energy 



internal 



h 2 nc , _ 

h hm > 

|k-ko|<K,< 



2U2 



h 2 k 
2m 



h 2 c 



2m(k - k o y 



(40) 



and C = limk.kQ^oolk — k | 4 n kcr . These two formulas are equivalent to Eqs. (j2J) 
and (j3J), in the contact-interaction limit. But now k can exceed 1/tq without 
leading to any problems, provided that the typical relative energy of these 
fermions is much smaller than H 2 /(mr^). 



5 Many-Body Problem in Momentum Space 



5.1 Basic Formulation 



Many-body theories are often studied in momentum space, because in the ther- 
modynamic limit, we often have translational symmetry, and the exploitation 
of this symmetry in momentum space simplifies many things. 

The s-wave contact interaction problem, however, lacks a satisfactory momen- 
tum representation to date. The peculiar contact interaction causes ultravio- 
let divergence problems, which some authors deal with by using concepts like 
"bare" coupling constants and renormalized ones; the shortcoming of these 
approaches is that the bare constants are ill-defined divergent quantities, and 
the sum of them and the divergent counterterms are ambiguous. 

Here we write down the momentum formulation of the problem, using the 
short-range selectors. 

We start from the coordinate space. The many-body Hamiltonian is 

H = internal + / dV £ V ext (r ) ^ (r ) fa (r ) , (41) 
internal = ~ / £ fa{v) V 2 fa{v) 

2m J a 

+ / d 3 rdVA(r')^(r)^{(r)^(r _ r 7 2 W r + r'/2). (42) 
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Note the similarity between the above equation and Eqs. (1) and (2) of 
Ref. [16j, where the hard-sphere Bose gas is studied; but in our context, in 
the contact interaction limit, the pseudopotential Hamiltonian is exact. Any 
right-hand side eigenstate |0) of the Hamiltonian (satisfying H\<f>) = E\(f))), or 
any linear combinations of such eigenstates (with an upper bound in energy), 
automatically satisfy the short range boundary condition 



dV^r'^jO - r'/2)^ T (r + r'/2)|0) = 0. (43) 

All the state vectors which satisfy this boundary condition form a subspace of 
the Hilbert space, and we call it the physical subspace, and represent it with 

v m. 



Equations (142]) and (143]) can be rewritten in the momentum space: 

-^internal = Y ^~7 C L C k CT + ~ q Y A ( k ') C q/2+kT C q/2-k| C q/2-k'| C q/2+k'T' ( 44 ) 
ka Z '" L qkk' 

^^( k ) c q/2~k|C q/2 +kTl0) = (for any q). (45) 

k 

All the ultraviolet divergence problems disappear, when A(k') is restored in 
the Hamiltonian. No divergent bare constants, no ad hoc regularizations, no 
renormalizations. Everything can be formulated simply and exactly. 

We now describe how to compute the expectation value of H internal under any 
state |0) G V. The rule is very simple: 

f k 2 

(^internal) = Y } Y 7717 ( C k CT C k<r) 



k 



2m 

Ana 



+ 



mQ qk , 



X^-M k ')( C q/2+kT C q/2-k| C q/2-k'| C q/2+k'T) f> ( 46 ) 



that is, the two terms should be grouped in the above way. Now the energy is 
finite, so the summation over k is convergent and the summand decays faster 
than 1/k 3 at large k. Consequently, we can insert the ^-selector, according to 
(1230): 

r k 2 

(^internal) = Y ^( k ) \Y TTZT ( C L C ka) 



2m 



qk' 

Distribute rj (k) to the two terms, we get 



+ _0 Y A ( k/ )( C q/2+kT C q/2-ki C q/2-k'| C q/2+k'T) 
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k 2 

(//internal) = &( k ) 7jH C L C kcr> 
kcr " 

C q/2-ki C q/2-k'| C q/2+k'T 

qkk' 

but the second term must vanish, since the Hermitian conjugate of Eq. (1451) 
is Ek^(k)(0| c q /2+kT c q/2-k| - °- Therefore 

(//internal) = E'iW^fcj- 

which is precisely the energy theorem. 

It is tempting to use this approach as the simplest way to derive the energy 
theorem. However it is not explained why the terms should be grouped in 
the way as in Eq. ( 1461) . even though it is a very natural grouping. It is the 
logical proof of the energy theorem in the last section that supports this way of 
grouping the terms. 

The inner product of Eq. ( 1451) with any quantum state is zero. This gives us 
many useful identities. The simplest of them is 

5>(k)(c_ ki c kT > = 0, (47) 

k 

which is a constraint on the pairing amplitude in the superfluid state. We 
have also seen the constraint on the two-body reduced density matrix, and 
we exploited it to rederive the energy theorem above. Similar constraints are 
present on the three-body, four-body, reduced density matrices. 



The internal Hamiltonian (1441) . although not Hermitian in the whole Hilbert 
space, is Hermitian in the physical subspace. This can be easily shown from the 
energy theorem, which indicates that the expectation value of //internal under 
any state in V is real. For any |0i), |0 2 ) e V, and any angle 9, |0i) + e l9 \(f>2) 
is also in V, and the expectation values of //internal under all these states are 

real. It then follows that (01 1 //internal 1 02 ) = (02 1 //internal) 01 )*• 

One might worry about the divergence of a in the unitarity limit, in which the 
second term on the right-hand side of ( 1441) appears ill-defined. We can show, 
however, that there is no real problem here. Before a becomes divergent, we 
should use the physical subspace condition to replace A(k') with A{k')+jr](k'), 
where j is an arbitrary constant. Choosing j = — 1 we get 

k 2 ^ 1 x| 

-f/internal|0) = ^ C kcr C kcr 1 4>) ~ ~ 7=7 - f/ ( k ') C q/2+kT C q/2-ki C q/2-k'| C q/2+k'T I <t>) > 



kCT 2m— mfi qkk/ 



(4£ 
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and the second term is now explicitly well-behaved in the unitarity limit. This 
equation is valid for any scattering length, no matter finite or divergent. To 
calculate the inner product of this equation with any state vector, we should 
use Eq. ( II 5 aft to convert the integral over k' to a limit at large k'. 

Note that in the unitarity limit the 77-selector is equal to the A-selector, making 
many expressions particularly simple. 



5.2 Few-body physics: an example 



All the few-body physics is clearly contained by our second-quantized formu- 
lation. 

To show that the formalism presented in this paper is a streamlined working 
method (and not only a formal framework), and to show that few-body physics 
and many-body physics can be treated with this very same formalism, we 
demonstrate a simple example, the low energy scattering of a fermionic dimer 
and a free fermion, in which case a > necessarily. This calculation can be 
compared with Ref. |25J in which the same problem is studied. 

It is desirable to work in momentum space because of spatial translational 
symmetry. The system is in a quantum state 

|0)= E 0k lk2k 34 lT i 2 T c k3|l O >> ( 49 ) 

kik 2 k 3 

where the wave function is antisymmetric under the exchange of ki and k2, 
and is zero if ki + k 2 + k 3 7^ 0. Substituting (H9|) into Schrodinger equation 
H |0) = E\(f)), where H is given by (j44l) (assuming V ext = 0), we get, 



1 ' 2 ' 3 j^j 



2m 
+ 



- > kik 2 k 3 



17lQ 

We define 



E A ( k i - K)<Pk^k + E A ( k ' 2 - k 3 )0k lk ^ 



k^kg kjkg 



0. (50) 



/k 2 = ^ E A ( k 'i - k 3 )0 k;k2k ,, (51) 



which is a function of just the length of k 2 , in the case of low energy s-wave 
scattering. If the relative momentum of the dimer and the fermion approaches 
0, E — —l/(ma 2 ) < 0, and Eq. (I5T)|) can be formally solved: 

, _ 47ra(/ kl - / k2 ) 

<? kl k 2 k 3 (jfe 2 + k 2 + fe 2)/ 2 + l/ a 2 ' W 
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for ki + k 2 + k 3 = 0. Substituting this back to the definition of /k 2 , and writing 
k 2 = k, k[ = -k/2 + k', and k^ = -k/2 - k', we get 

/k = to /g A(k0 _^_^ , (88) 

where the integral of the second term on the right-hand side can be immedi- 
ately evaluated, using fl9e|) . In the first term, however, A(k') can be dropped, 
because the integrand decays faster than 1/k' 3 for k' — > oo. If this were not the 
case, /k would have to decay like 1/k or even slower at large k, and according 
to Eq. (J52D, 0k! k 2 k 3 would have to decay like 1/K S or slower, when ki, &2 and 
k% are all of the order K — > oo, making the wave function unnormalizable at 
short distances [26]. We thus get 

i ^ 7T ^- a ' 1 )/^^/|| ^ +k . k f k ; <:2+a - 2 -°. («) 

where k' has been shifted. 



The Fourier transform of /k is proportional to the relative wave function of the 
fermion and the dimer at large distances, and must be of the form 1 — cifd/r 
for r — > oo. Here a/d is the fermion-dimer scattering length, whose value we 
shall determine. Consequently, /k is of the form (27r) 3 #(k) — Anafd/k 2 plus 
regular terms at small k. 

Let /k = F(t), where t = (ka) 2 is dimensionless. Equation ( 1541) is easily 
simplified to 

(Jl + — -l) F(t) + -^7= r Fit 1 ) arctanh ^ dt 7 = 0. (55) 

It can be shown that F(t) oc —irS(t)/\/i + (a^/a)/t plus a smooth function 
for all t > (including the neighborhood of t = 0), and that F(t) oc cot^/ 2 + 
Ci £(mo-i)/2 pi us hig^gj; order terms at t — > 00, where 



= -4.16622197664779257337, and c x /c = 0.30268913080233667524. 

Exploiting these properties, we can discretize the integral equation in an ap- 
propriate way, and solve it extremely accurately. The resultant fermion-dimer 
scattering length turns out to be 



a fd = 1.1790662349 a. (56) 

To make an extremely accurate prediction for a real system (eg, ultracold 
atoms), we have to take various corrections to the idealized s-wave contact 
interaction into account, a topic beyond the scope of this paper. 
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6 Summary and Outlook 



To solve a paradox arising in the formulation of the s-wave contact inter- 
action problem of ultracold fermions in momentum space, we modified the 
conventional notion of integrals. This unexpectedly led us to a simple relation 
between the energy and the momentum distribution of these fermions. 

The new notion is that while the integral of an ordinary function over the 
whole space is always equal to the limit of integrals over finite regions of 
space, some generalized functions do not obey this rule. For them, certain 
full-space integrals are specified beforehand, in a way that is compatible with 
their other properties and with conventional mathematics. 

Our approach is very different from the ones used by some people in some rela- 
tivistic quantum field theories (in which they force integrals like / d A kk~ 2 = 0), 
even though one might form this false impression if one only reads Sec. |2] of 
this paper. In our approach, for example, A(k)/k 2 and r](k)/k 2 are equal for 
all finite k's, but their full-space integrals are different. The infinitely many 
possibilities form a linear space, called the selector space, in which each point 
corresponds to certain values of some integrals. Also, in our approach, inte- 
grals of ordinary functions which are divergent remain divergent forever (for 
example, / d 4 kk~ 2 = oo), and the generalized functions and ordinary ones 
coexist, forming a coherent system. 

We discussed the physical implications of our calculations. In particular, if we 
measure the momentum distribution of fermions with s-wave contact interac- 
tions, and plot the result in a certain way, a simple asymptote appears, and 
the height of such an asymptote at a certain horizontal coordinate is equal 
to the internal energy. This result is the generalization of a simple property 
of noninteracting particles, in the context of strongly interacting ultracold 
quantum gases with resonant interactions (\a\ ^> r ). 

We also discussed the application of the results to confined quantum gases, 
and found that the energy of the system can be extracted from the equal-time 
one-body reduced density matrix only. This is of course completely different 
from the density functional theory of electrons, in which the explicit form of 
the functional is unknown. 

We formulated the s-wave contact interaction problem of ultracold fermions in 
a simple and coherent way, and can now freely transform the problem between 
coordinate space and momentum space. This is a streamlined formalism for 
studying many concrete problems, including few-body ones. 

Full applications of this formalism to the many-body problem are a topic 
of future research. Because of the peculiar structure of this system, we can 
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not predict whether there will be further surprises awaiting us. The short- 
range boundary condition of the many-body wave function, essentially the 
definition of the scattering length, appears to be a highly nontrivial constraint. 
Can we construct a series of approximations to reduce discrepancies with this 
constraint in a progressive way, and to minimize the system's free energy at 
the same time? It is possible. 

The BCS wave function is of course not consistent with this constraint. For 
example, Ek' ^( k ')(0BCs|Cq /2+kT c] l/2 _ ki c q/ 2_ k ^c q/2+k / T |0BCs) ^ at q ^ 0. 

I hope to extend some essential ideas in this paper to all systems with contact 
interactions, including the Standard Model of particle physics. If there are 
R independent "renormalized" constants that must be determined by exper- 
iments (or by a higher-energy theory) in a low-energy effective theory, is the 
dimension of the associated selector space equal to R + 1? We have this tenta- 
tive guess, because in the nonrelativistic s-wave contact interaction problem 
considered in this paper, we have only one "renormalized" constant, the scat- 
tering length, and our selector space is two-dimensional. Another conjecture 
is that the selector spaces of those other theories may have some nontrivial 
algebraic properties; for example, the momentum translation of one selector 
may lead to mixtures with other ones. It may also be interesting to investi- 
gate how these selectors evolve (inside the selector space) under the continuous 
scaling transformation, to see whether or how the conventional notion of renor- 
malization group flow is incorporated. In the two-dimensional selector space 
described in this paper, A(k) is unchanged if k is rescaled, but L(k) undergoes 
a simple scaling transformation, and a generic selector [linear combination of 
A(k) and L(k)] generally changes after a scaling transformation. 

If such a universal approach is realized, certain developments will follow. The 
most important of them will probably be related to some nonperturbative 
properties of interacting quantum fields. Even if we can not determine all 
the properties quantitatively, we may still find some exact relations and/or 
qualitative features. 
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